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1. Introduction 

Determination of the leading-order parameters of chiral perturbation theory (ChPT), i.e. the 
chiral condensate £ and the pion decay constant F, is the first goal for lattice QCD in establishing 
the link between QCD and ChPT. This is not an easy task, because the large volume limit must be 
taken before taking the chiral limit. Both limits require enormous computational cost. 

Recently, another approach is getting popular, that is the lattice calculation in the £-regime. 
Here, one reduces the quark mass very close to the chiral limit at a fixed volume V. Of particular 
interest is the region where the pion correlation length, or the inverse pion mass 1 jm % , overshoots 
the size of the box L. In this regime, the finite volume effect becomes prominent, but it can be 
treated in a systematic way within ChPT as an expansion in a new parameter e [Jl|, ||, [§], which 
has the scale m n /K cui ~ p 1 / 'A^ ut ~ £ 2 , where p denotes the pion momentum and A cut is the cutoff 
of ChPT. Zero-momentum mode of pion is treated non-perturbatively at the leading order of the 
£-expansion. With this expansion, precise predictions for the volume and topological charge de- 
pendences of the quark condensate, meson correlators, etc. are obtained in terms of the low-energy 
constants in the infinite volume. 

Simulating lattice QCD in the £ -regime has potential advantages. Computational cost can be 
kept manageable by staying in a small box, when reducing the quark mass until the chiral extrap- 
olation becomes unnecessary. The remarkable topological charge and quark mass dependences of 
the meson correlation functions, for instance, in the £ -regime are helpful to identify the correspon- 
dence between QCD and ChPT. These nice properties have already been exploited in the quenched 
QCD studies [Q]. In such works the overlap-Dirac operator ||5|] is used, since any tiny violation of 
chiral symmetry may be amplified in the £-regime. 

A recent work of the JLQCD collaboration [§] extended the study in the £-regime to un- 
quenched QCD. (See also [Q] for smaller scale works.) On a 16 3 x 32 lattice at the lattice spacing 
a ~ 0.1 1 fm (determined from ro ~ 0.49 fm) we have generated 460 configurations with two-flavor 
dynamical overlap quarks and Iwasaki gauge action. We carried out a simulation at an extremely 
light sea quark mass m ~ 3 MeV, which is within the £-regime. Comparing the Dirac spectrum with 
the predictions of chiral Random Matrix Theory, we extracted the value of £ at the leading order 
of the £-expansion as £ MS (2GeV) = [25 1(7) (11) MeV] 3 , where the second error is an estimate of 
the systematic error due to the NLO effects of the £-expansion. 

In this work, we use the same configurations in the £-regime to calculate the meson correlators 
in various channels. Analytic predictions of ChPT for various channels are known to the next-to- 
leading order of the £-expansion [Q, [Kj. Their extensions to the partially quenched ChPT have 
also become available recently [111]]. We compare these ChPT predictions with the lattice data to 
extract L and F to the NLO accuracy. 



2. (Partially quenched) chiral perturbation theory in the e-regime at fixed topology 

First, we briefly review the results for the meson correlators calculated within (partially quen 
-ched) chiral perturbation theory. See [|9|, [K], 11] for details. Here we consider 7V V valence quarks 
with a mass m v and Nf = 2 sea quarks with a mass m s , both in the £-regime. 
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As a fundamental building block of this section, let us define the partially quenched zero-mode 
partition function which consists of two physical quarks (of which mass is m s ) and one valence 
quark (with m v ) and 1 bosonic quark (with nib), at a fixed topology v. 



^v pq (^IMv,m, 



(A 2 -A 2 ) 2 



det 



/ Ky(n b ) I v (n v ) I v (n s ) Iy^Hs)/^ 
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where = m^Ey, /x v = m v LV and /i v = m s £V. ^ v 's and 7 v 's are the modified Bessel functions. 
Note that in the limit ju& — ► /x v , it reproduces the Nf = 2 full theory partition function 



^VvlMv,^) = ^(lO - det f 



r) 



(2.2) 



The partially quenched chiral condensate at finite V and v is defined by 



= - hm 3 — hi&v^iUbl^Hs). 



(2.3) 



It is not difficult to see that in the jx v — > jx s limit, the partially quenched condensate becomes the one 
in the full theory, = &hl = 1 * injf* 11 ^). In t h e following, we also use the double 

derivative of the condensate, defined by 



AEv Q (Mv,Ms) _ lim M ^ v d^,jf Q (Mfo|Mv,jU.v) 



(2.4) 



First, we present the correlation functions of the flavored pseudo-scalar and scalar operators, 
P a (x) = q(x)T a y 5 q(x) and S a {x) = q(x)z a q(x), where z a denotes the generator of SU(N V ) group 
which the valence quark field q(x) belongs to. For these correlators, the partially quenched expres- 
sion at fixed topology is known to &{e 2 ) (no sum over a) [11 1: 



C P {t) = d 3 x(P"(x)P"(0)) 
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C s (f) = / J 3 x(S fl (x)S"(0)) 
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h{t/T), (2.6) 



where jU? ff = ra,£ e ftV, and the ? dependence of correlators is represented by a function h\{t/T) = 
(t/T — 1 /2) 2 — 1/12 /2. £ receives one-loop correction as £ e ff = E (l + 3p\/2F 2 V 1 ^ 2 ) , where 
j8i is the so-called shape coefficient. In our numerical study, j8i = 0.0836. 
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Next, consider the flavored axial-vector and vector operators Aq(jc) = q{x)T a yoysq{x) and 
V$(x) = q(x)l a Yoq{x). Their correlators at <^(e 2 ) in the N f = 2 ChPT are (see |[k]] for details) 



2 / J81 


n 






4^4 ull (Ms) 


T 2 


F 2 L 


V 


2 / J8i 


go 
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+ «_v_^ —h{t/T) \ , (2.7) 



(2.8) 



where &oo is another numerical factor depending on the shape of the box (in our case, &oo 
0.08331). are defined by 



f0- 



= > 3 T 1±2 



ft* £ (udt)2 



,(2.9) 



Note that C P (t) and C s (f) are sensitive to I and F enters only at &(z 2 ). For C A (0 and C v (t), 
on the other hand, F appears in the leading term and thus can be extracted efficiently. 

3. Lattice simulations 

Here we summarize our numerical set up for the simulations. For other details, see Our 
lattice size is 16 3 x 32 and the lattice spacing is determined as a = 0.1111(24) fm from the heavy 
quark potential assuming ro = 0.49 fm. We use the overlap fermion of which the Dirac operator 
with a quark mass m is given by 

m + -j + [mo - 2) »sgn[ff w (-mo)], (3.1) 

where Hy?{— mo) = Y5D w (—m ) denotes the standard Hermitian Wilson-Dirac operator. We choose 
mo = 1.6 throughout this work. (Here and in the following the mass parameters are given in the 
lattice unit, unless otherwise stated.) For the gauge action, we use the Iwasaki action at j8 = 2.35 
together with an additional determinant factor corresponding to Wilson fermions and associated 
twisted-mass ghosts [Q], which forbids the topology changes along the Monte Carlo updates. In 
this work, topological charge is fixed to v = 0. We use the hybrid Monte Carlo (HMC) algorithm. 



The sign function in (3.1 ) is approximated by a rational function with Zolotarev's coefficients after 
projecting out a few lowest-lying eigenmodes. 

In this work, we focus on the run at the lightest sea quark mass m = 0.002, which corresponds 
to ~ 3 MeV in the physical unit and the system is well within the £-regime. We accumulated 
4,600 trajectories after discarding 400 trajectories for thermalization. At every 10 trajectories, we 
calculated the meson correlators in various channels. We take four values of valence quark mass, 
m = 0.0005, 0.001, 0.002, and 0.003 (1-4 MeV). The inversion of the Dirac operator is performed 
simultaneously for all the valence quark masses using the multimass solver. 

For the meson correlators, the purely low-mode contribution (both quark and anti-quark propa- 
gators are represented by the low-lying eigenmodes) is replaced by the low-mode averaged (LMA) 
one, i.e. the source point is averaged over all lattice sites. We find that LMA with 100 eigenmodes 
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improves the statistical signal substantially for the pseudo-scalar and scalar correlators while the 
improvement is marginal for the axial-vector and vector correlators. 

The jackknife bin-size is chosen as 20, with which the statistical error saturates. 



4. Numerical results 



Now we present the numerical results. 

For the axial current we use the local operator Aq(x), which receives finite renormalization. 
We calculated the renormalization factor using the RI/MOM scheme and obtained = 1.3513(13). 
In the following, we neglect this tiny statistical error for Z^. 

Fig. |] (left panel) shows the data for the axial-vector correlator at m v = m s = 0.002. A 
two-parameter fit to (2/7) works well (# 2 /d.o.f ~ 0.01) for the low-mode averaged data (filled 
squares) as shown by a blue curve. The fitting range is t G [12,20]. From this fit we obtain 
£ ~ [260(32) MeV] 3 and F ~ 90(6) MeV. The statistical error in £ is large (~ 30%) because £ 
appears only at 6{z 2 \ while F is determined to a good precision. 

Next, let us look at the pseudo-scalar channel at m v = m s = 0.002. Using F obtained via 
the axial-vector correlator as an input, we fit the pseudo-scalar channel with ( |2~5| ) and obtain £ 
with much better precision. Fig. p] (right panel) shows the fit curve with a fit range t £ [12,20] 
(# 2 /d.o.f = 0.07). A two-parameter fit with F as another free parameter does not work because the 
sensitivity to F is too weak. 

To obtain the best result we perform a simultaneous fit to both P a P a and A[|Aq correlators and 
obtain I = [227.6(3.7) MeV] 3 and F = 87.3(5.6) MeV (or y/lF = 123.5(7.9) MeV), where the 
statistical error in a is also taken into account. Here the fit range is t G [12,20] and # 2 /d.o.f = 0.02. 
Multiplying the renormalization factor calculated non-perturbatively through the RI/MOM scheme, 
we obtained the renormalized condensate £ MS (2 GeV) = [239.8(4.0) MeV] 3 . This value is consis- 
tent with our previous result obtained through the Dirac spectrum £ Dirac = [25 1 (7 ) ( 1 1 ) MeV] 3 . 

Once the parameters £ and F are determined, there is no additional free parameters at the 
given order of the £-expansion, and therefore the comparison provides a stringent test of the lattice 
data and/or the £-expansion. For instance, the lattice data for the scalar channel is shown in Fig. || 
(left panel). The curve in the plot is not a fit to the data but a ChPT prediction (|l6|). The agreement 
is remarkable. Furthermore, we also test the consistency with the partially quenched data sets 
(m v ^ m s ) for the pseudo-scalar channel (Fig. || (right panel)). The curves showing the ChPT 
prediction without free parameters are perfectly consistent with the lattice data. 



5. Discussions 



We have compared the different channels (pseudo-scalar, scalar, axial-vector and vector), par- 
tially quenched correlators with four different valence quark masses, and Dirac spectrum, and they 
are all consistent with each other. For the change of the fitting range within ? m i n € [10, 15], both 
£ and F are quite stable (within 1%) with similar error-bars. For the finite volume correction, the 
NLO in the e-expansion, i.e. 6{e 2 ) ~ p 2 ~ 1 /L 2 , is taken into account within ChPT. 

Still, we find some discrepancy in the pion decay constant with the preliminary result (F = 
78(3)(1) MeV) obtained in the /^-regime using the NNLO formula for the chiral extrapolation [|l"2|]. 



5 



Meson correlators in the £-regime 



H. Fukaya 



CM 



0.0001 



8e-05 



6e-05 



4e-05 



2e-05 



with LMA ma=0.002 
withoutLMA ma=0.002 
ChPT (£=[260MeV] 3 , F =89.5MeV.) 




10 15 20 25 30 
t 



0.01 



lattice ma=0.002 ■ — - 
ChPT (Z=[227MeV] , F,=89.5MeV.) 

0.008 I I 



0.006 



0.004 



0.002 




1 1 1 1 1 1 ■— 1 

5 10 15 20 25 30 

t 

Figure 1: Upper figure: the axial vector correlator in the £-regime. Filled square symbols denote the low- 
mode averaged correlator while the crosses are not averaged. The solid curve shows the best fit yielding 
£ ~ [260 MeV] 3 and F ~ 90 MeV). Lower figure: the pseudo-scalar correlator. The solid curve shows the 
best fit yielding E - [227 MeV] 3 (F = 89.5 MeV is given as an input). 



(The /^-regime simulation has been done at a slightly coarser lattice spacing, a ~ 0.12 fm, but we 
do not expect substantial discretization effect by this small change of a.) This may signal some un- 
known source of systematic error, possibly the higher order contributions in either the £-expansion 
or the /^-expansion. 

HF thanks P. H. Damgaard, T. DeGrand and P. Hasenfratz for fruitful discussions. Numeri- 
cal simulations are performed on Hitachi SRI 1000 and IBM System Blue Gene Solution at High 
Energy Accelerator Research Organization (KEK) under a support of its Large Scale Simulation 
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Figure 2: The scalar (left) and the partially quenched pseudoscalar (right) correlators in the e-regime. The 
solid curves represent the ChPT prediction with £ = [227.6MeV] 3 and F = 87.3MeV (No free parameter 
left). 
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